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PREDICTION  OF  SOLDER  JOINT  FATIGUE  LIFE 


H.D.  Solomon,  V.  Brzozowski,  and  D.G.  Thompson 


1.  INTRODUCTION 

The  phenomenon  of  the  fatigue  of  surface  mount  joints  has  been  the  subject  of 
considerable  study  in  recent  years  [1-11].  This  present  study  aims  at  predicting  the 
isothermal  low-cycle  fatigue  (LCF)  of  a  leadless  Chip  Carrier/Printed  Wiring  Board 
(CC/PWB)  joint  from  isothermal  LCF  data  and  a  finite  element  analysis  of  the  strains 
developed  in  the  joint.  This  prediction  is  compared  with  fatigue  lives  measured  in 
tests  on  actual  leadless  CC/PWB  joints.  These  joints  were  made  using  typical 
manufacturing  procedures  and  actual  chip  carriers  and  printed  wiring  boards. 

The  procedures  utilized  in  the  calculation  are  discussed  in  detail  as  are  any  approx¬ 
imations  utilized  and  the  methods  used  to  correlate  the  prediction  with  the  experimen¬ 
tally  determined  fatigue  life.  Of  particular  interest  is  a  size  effect  which  must  be  con¬ 
sidered  when  applying  the  LCF  data  to  actual  joints.  This  size  effect  results  from  the 
behavior  of  thin  solder  layers  tested  in  shear  as  compared  to  more  typical  tests  on  bulk 
specimens  tested  in  tension.  The  most  important  aspect  of  this  paper  is  not  the 
specific  correlations  obtained  between  measured  and  calculated  fatigue  lives.  Rather, 
what  is  most  significant  is  the  description  of  the  approach  that  was  followed,  as  this 
approach  enables  one  to  apply  the  general  isothermal  fatigue  life  to  joints  with 
different  geometries.  Thus,  this  approach  is  described  in  considerable  detail. 

This  work  was  done  as  part  of  a  U.S.  Air  Force- sponsored  manufacturing  technol¬ 
ogy  program,  with  participation  by  Martin  Marietta  Orlando  Aerospace,  Westing- 
house  (Defense  Electronics  and  R&D  Center)  and  GE  (Avionics  and  Electronics  Sys¬ 
tems  Division  and  R&D  Center).  The  finite  element  analysis  of  solder  joint  strains 
was  performed  by  Westinghouse  personnel.  The  fatigue  tests  on  CC/PWB  joints, 
analysis  of  the  LCF  data,  development  of  techniques  to  predict  the  fatigue  life,  and 
correlation  of  the  experimental  data  with  life  predictions  was  done  at  the  GE  R&D 
Center. 

1.1  Procedures 

No  new  experimental  data  is  introduced  here,  so  the  reader  will  be  directed  to  the 
data  sources  for  the  specific  details  utilized  in  the  tests  to  acquire  that  data.  It  is  im¬ 
portant,  however,  to  get  an  overview  of  what  is  being  done,  as  several  different  data 
sources  are  utilized. 

The  calculation  of  the  fatigue  life  is  divided  into  two  distinct  parts:  the  calculation 
of  the  strains  developed  in  the  joints  and  the  prediction  of  the  fatigue  life,  Nf,  from 


these  strains.  The  strains  were  determined  from  a  finite  element  analysis  of  solder 
joints  [12].  The  fatigue  life  -  strain  correlation  was  made  from  isothermal  tests  on 
simple  solder  joints  [2-5]  (but  not  CC/PWB  joints).  This  predicted  fatigue  life  is  com¬ 
pared  to  that  measured  on  actual  CC/PWB  joints  [13].  The  FEA  was  done  on  joints 
which  had  the  same  geometry  as  the  tested  joints  or  bracketed  the  critical  feature  of 
the  amount  of  solder  in  the  fillet. 

A  simple  flow  chart  illustrating  the  approach  taken  here  is  shown  in  Figure  1.  The 
FEA  uses  a  displacement  A  and  calculates  the  strain  distribution  -77.  These  strains  are 
in  turn  used  to  calculate  Ay  via  a  Coffin-Manson  LCF  curve  developed  from  the  tests 
on  simple  solder  joints.  This  life  is  compared  to  tests  run  [13]  on  actual  CC/PWB 
joints  which  were  subjected  to  the  same  sort  of  displacement.  A,  as  was  used  in  the 
FEA. 

The  analysis  is  complicated  by  the  fact  that  the  actual  strain  distribution  in  a  solder 
joint  is  highly  nonuniform.  Without  a  simple  77,  it  is  necessary  to  use  a  cumulative 
damage  sort  of  approach  to  calculate  A y  for  various  conditions.  The  use  of  the  LCF 
curves  is  further  complicated  by  a  size  effect  which  must  be  corrected  for  before  it  can 
be  used  for  the  prediction  of  actual  solder  joint  lives.  The  last  set  of  complications 
comes  in  utilizing  the  proper  displacements  when  correlating  A  for  the  FEA  and  for 
the  direct  measurement  of  the  joint  fatigue  life.  The  FEA  utilized  a  zero  to  A  dis¬ 
placement,  whereas  the  fatigue  tests  were  run  from  -A  to  +  A.  Furthermore,  the  A 
utilized  in  the  FEA  was  a  total  displacement,  whereas  plastic  displacements  were  uti¬ 
lized  in  measuring  Ay  for  the  solder  joints.  It  is  therefore  necessary  to  correct  for  the 
elastic  strains  and  this  will  be  discussed. 

I.  Finite  Element  Analysis 

The  FEA  model  used  for  analyzing  leadless  chip  carrier  solder  joint  stresses  and 
strains  induced  by  thermal  or  mechanical  excursions  was  developed  by  Westinghouse 
under  the  Air  Force  contract  F-33615-82-C-5047  supported  by  the  Wright  Patterson 
Materials  Laboratory. 

A  detailed  solder  joint  FEA  model  allows  the  user  to  define,  in  great  detail,  the 
outer  dimensions  of  one  solder  joint  and  then  generates  a  finite  element  mesh  of  the 
solder  volume.  The  model  assumes  that  the  solder  is  attached  to  the  HCC  and  that 
this  interface  is  a  fixed  boundary.  A  copper  pad  is  attached  to  the  bottom  surface  of 
the  solder  and  any  deflection  due  to  mechanical  or  thermal  loads  is  applied  to  the 
copper  pad.  The  mesh,  material  properties  of  the  solder  and  copper,  boundary  condi¬ 
tions  and  loading  conditions  are  input  to  WECAN  [  17]. 

WECAN  calculates  the  detailed  strain  distribution  in  the  solder  joint  in  question 
under  a  set  of  given  loading  conditions  and  allows  the  user  to  study  the  levels  of  strain 
and  their  location. 
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FEA  LCF 


In  N, 

N<  measured  as  a  function  of  A 


Figure  1.  Approach  to  relate  joint  strains  determined  by  FEA  with  strain  life  LCF  data 

and  predict  the  joint  fatigue  life. 


Visual  x-ray  and  cross-sectional  inspection  techniques  were  used  to  characterize 
the  solder  joint  geometries  fabricated  for  the  LCF  tests.  Due  to  manufacturing  pro¬ 
duction  variations  incurred  during  sample  fabrication,  solder  joint  geometries  varied 
from  sample  to  sample  as  well  as  from  joint  to  joint  within  a  sample.  The  wider  varia¬ 
tions  in  solder  joint  geometry,  due  to  these  production  variations,  make  it  impossible 
to  analyze  every  possible  configuration.  The  decision  was  made  that  for  this  set  of 
analysis,  the  modeling  would  concentrate  on  two  factors,  solder  joint  height  and  fillet 
shape.  For  the  solder  joint  height  (height  between  the  top  of  the  copper  pad  on  the 
PWB  and  the  metallization  on  the  HCC),  heights  of  2  mils,  3  mils,  and  5  mils  were 
chosen  to  bracket  the  cross-sectional  heights  reported  in  the  microsection  work  with 
the  3  mil  height  acting  as  a  nominal  or  baseline  height.  Two  solder  joint  fillet  sizes 
were  chosen,  a  bulbous  or  full  fillet  and  a  minimal  fillet,  shown  in  Figures  2  and  3. 

The  elastic  properties  required  for  the  copper  and  solder  used  in  the  analysis  are 
contained  in  Table  1.  The  plastic  properties  of  the  solder  (i.e.,  stress/strain  curves) 
used  are  given  in  Figure  4.  These  curves  were  developed  from  bulk  63-37  solder  speci¬ 
men  tensile  tests  performed  at  Westinghouse  under  contract  number  F-33615-82-C- 
5047.  The  testing  involved  a  cast  ingot  of  solder  which  underwent  isothermal  cyclic 
testing  [16]. 

The  FEA  used  data  obtained  from  tests  on  actual  Chip  Carrier/PWB  joints  to  set 
the  parameters  of  the  model.  These  were  the  same  specimens  used  for  the  isother¬ 
mal  LCF  tests  [13].  The  loading  of  the  specimens  was  nominally  in  simple  shear  and 
this  was  the  assumed  loading  for  the  FEA.  The  total  elastic  plus  plastic  deflection 
measured  by  the  LCF  test  extensometers  within  the  first  few  cycles  of  the  test  was 
used  as  a  basis  to  calculate  the  impose  deflections  used  in  the  analyses. 

It  was  assumed  that  the  measured  deflection  was  equally  divided  between  each 
row  of  solder  joints  and  that  the  deflection  was  equal  in  both  the  tensile  and  compres¬ 
sive  load  cycle,  ( P  +  and  P').  The  imposed  deflection  in  each  direction  on  any  one 
solder  joint  during  the  analysis  was  thus  calculated  by  taking  1/4  of  the  total  deflection 
from  tip  to  tip  of  the  elastic  and  plastic  hysteresis  loop.  In  addition,  each  solder  joint 
model  analysis  had  attached  to  the  model,  in  the  direction  of  the  loading,  a  spring  with 
a  very  large  spring  constant  (see  Figure  5A).  The  deflection  imposed  on  the  solder 
joint  causes  the  spring  to  deflect  and  the  deflection  of  the  spring  times  the  spring  con¬ 
stant  gives  the  loading  on  the  spring.  The  high  spring  constant  was  chosen  such  that 
the  deflection  in  the  spring  was  much  less  than  the  deflection  in  the  solder  joint.  Since 
the  spring  is  in  parallel  with  the  solder  joint,  this  is  the  loading  applied  to  the  solder 
joint  which  deflects  it  a  prescribed  amount. 
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Figure  3,  FEA  determined  strains  Tor  a  minimal  joint. 
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Table  1 

Copper  and  Solder  Elastic  Properties  Used  in  Computer  Analyses 


Copper 

Solder 

SN63 

Coefficient  of  Thermal  Expansion 

10.5  x  KT6 

15.7  x  10'7 

Poisson’s  Ratio 

.36 

.35 

Elastic  Modules 

17  x  106 

3.6  x  106 

Note: 

Since  the  HCC  to  solder  joint  interface  is  considered  to  be  a 
fixed  surface,  the  value  of  the  CTE  for  the  HCC  6  x  10* 
in/in/°C  must  be  subtracted  from  the  values  for  solder  and 
copper  when  performing  any  thermal  stress  analysis. 
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Figure  4.  Stress/strain  input  data. 
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Spring  Element  Determines 
Load  Solder  vs  Deflection  Solder 


^SP  >  >  ^  Solder 

Deflection  Spring  <  <  Deflection  Solder 
Load  Spring  =  Load  Solder 

(a) 


Figure  5.  a.  Spring  attachment  to  solder  joint,  b.  Correlation  of  LCF  test  data  to  modeling 
prediction. 
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Figure  5B  plots  the  computer  calculated  change  in  force  A P  vs.  total  deflection 
A(A+ B)  across  two  rows  of  11  solder  joints  which  are  3  mils  high  with  a  bulbous  fillet 
and  3  mils  high  with  a  minimal  fillet.  Superimposed  on  these  values  are  data  points 
from  samples  manufactured  by  GE  and  Westinghouse  which  had  undergone  the  LCF 
test.  The  calculated  force  values  bracket  the  actual  test  data  very  well.  The  data 
correlation  gives  confirmation  as  to  the  validity  of  the  FEA  model  and  the  calculated 
strain  in  the  FEA  solder  joint. 

Examples  of  the  results  of  the  FEA  analysis  are  given  in  Figures  2  and  3.  These 
1/2  models  detail  the  Von  Mises  strain  levels  induced  in  a  3  mil,  high  minimal  and 
bulbous  solder  joint  undergoing  a  total  deflection  of  .25  mils  in  the  y  direction.  An  in- 
depth  analysis  of  the  findings  of  the  FEA  Analysis  are  contained  in  Ref.  12.  The  fol¬ 
lowing  sections  of  this  paper  discuss  the  LCF  tests  and  correlate  the  measured  fatigue 
lives  to  the  induced  strain  levels  predicted  by  the  FEA  analysis. 

II.  Low-Cycle  Fatigue  Analysis 

The  calculation  of  the  total  shear  strain  7,  enables  one  to  estimate  the  fatigue  life. 
A  complication  arises,  however,  because  there  is  no  single  7,  which  defines  the  joint 
strain.  Tests  run  on  actual  joints  have  shown  that  the  initial  cracking  occurs  under  the 
chip  carrier,  where  the  strain  is  the  highest.  The  cracks  then  propagate  into  the  fillet 
where  the  strains  are  lower.  Figure  6  illustrates  this  schematically.  Cracking  under 
the  chip  carrier  occurs  due  to  the  application  of  strain  71 .  N 1  cycles  are  required  to 
produce  a  crack  length  L  j.  The  strain  in  the  transition  into  the  fillet  is  72,  requiring 
N2  cycles  to  run  a  crack  length  L2.  The  final  stage  in  failure  consists  of  growing  the 
crack  through  the  fillet  where  the  strain  is  73,  length  L3,  requiring  N 3  cycles. 

Table  2  lists  71 , 72  and  73  for  the  thin  and  bulbus  fillet  joints  analyzed  in  Figures  2 
and  3.  These  strains  are  the  estimated  weighted  averages  over  lengths  L\,L2  and  L  ?. 
These  strains  cannot  be  used  directly  to  estimate  N N2,  or  AI3  because  they  are  the 
VonMises  strains,  i.e., 

<tVM  =  -£2)2  +  (ft  '  £3 )*  +  fo  -f3)2)  1  PI 

where  *s  the  equivalent  VonMises  strain  and  e(,  z2  and  £3  are  the  principal  strains. 

The  fatigue  life  data  we  will  use  was  generated  in  shear  so  we  must  determine  a  shear 
equivalent  to  the  VonMises  strain.  For  pure  shear 

l~£\-^2  [-] 

C\  =  -£2.  £3  =  o  [3] 


* 
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CC/PWB  Joint  Failure  Process 


N,  Cycles  under  7,  strain  for  L,  distance 
N 2 Cycles  under  7,  strain  for  L2  distance 
N,  Cycles  under  7}  strain  for  L3  distance 
N,~25%  load  drop 
N,  +  N2  ~  50%  load  drop 
N,  +  Nj+  N  j~100%  load  drop 

Figure  6.  Schematic  representation  of  cracking  in  CC/PWB  joints. 


Table  2 

Strains  Developed  in  Joints 


(Under  cc) 

72  (Transition) 

73  (Fillet) 

Bulbus  (VonMises) 

0.050 

0.025 

0.0050 

Minimal  (VonMises) 

0.055 

0.030 

0.011 

Mean  (VonMises) 

0.0525 

0.0275 

0.008 

It  (Equivalent  Shear) 

0.091 

0.0477 

0.0139 

A-?;- 

0.182 

0.095 

0.028 
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1  =  2^!  =  -2e2 


[4] 


so,  from  Equations  1-4,  we  get 

£VM  =  7f 1  =  °'577'r  [5] 

The  vr  (equivalent  shear  strain)  was  calculated  by  dividing  the  VonMises  strains 
by  0.577  (as  per  Equation  5).  The  subscript  T  is  a  reminder  that  the  calculation  was 
for  the  total  (elastic  +  plastic)  strain.  The  mean  of  the  strains  in  the  Bulbus  and 
Nominal  geometries  is  being  used  because  the  actual  joints  were  between  these  two 
extremes.  As  can  be  seen,  there  is  only  a  modest  difference  in  the  strain  distribution, 
so  using  a  mean  value  is  very  reasonable. 

The  last  row  of  Table  2  lists  A77  which  is  just  2  77.  The  FEA  model  utilized  O  to 
A  displacements,  whereas  in  testing  actual  joints  a  -A  to  +A  displacement  range  was 
applied.  Symmetry  is  being  assumed  because  the  measured  [  13]  hysteresis  loops  for 
tests  on  actual  joints  are  symmetric,  so  A77  =  2  77- . 

The  LCF  data  utilized  the  plastic  strain  range  A7 P  to  correlate  the  fatigue  life  [2- 
5],  so  it  is  necessary  to  correct  this  data  to  yield  a  A77  -  Nf  curve,  but  first  we  need  to 
consider  the  A7p  -  Nf  curves  shown  in  Figure  7.  The  data  on  Figure  7  comes  from  a 
previous  study  [2]  of  single  solder  joints  tested  in  simple  shear  (simple  shear  is,  except 
for  a  rotation,  the  same  as  pure  shear,  and  Equation  (5)  still  holds).  In  this  study  a 
single  ~  0.0075  in.  (0.19  mm)  solder  layer  was  used  to  join  Cu  or  brass  blocks.  The 
area  soldered  was  0.1  in.  (2.54  mm)  by  0.5  in.  (12.7  mm)  with  the  shear  in  the  0.5  in. 
direction.  No  appreciable  fillets  were  made  with  these  joints  so  the  strain  could  easily 
be  calculated  by  dividing  the  applied  displacement  by  the  thickness  of  the  solder  layer. 

Two  curve  fit  correlations  are  shown  in  Figure  7.  The  correlation  with  a  =  0.52 
(the  Coffin- Manson  LCF  exponent)  has  been  discussed  previously  [2].  The  slope  was 
determined  from  a  least  squares  fit  with  a  conventional  statistical  computing  package 
which  assumes  the  ordinate  (A7 p)  to  be  the  dependent  variable  and  curve  fits  to  mini¬ 
mize  the  deviations  from  the  ordinate  coordinates.  Unfortunately,  the  convention 
with  fatigue  data  is  to  plot  the  independent  variable  A7 p  on  the  ordinate.  The  solid 
line  (with  a  =  0.59)  was  calculated  by  correctly  minimizing  the  deviations  from  the 
dependent  Nf  values.  Over  the  range  of  the  data,  there  is  little  difference  between  the 
two  curves  and  both  give  essentially  the  same  life  predictions.  Extrapolating  the  data 
to  lower  values  of  A7 p  does,  however,  lead  to  differences.  Differences  also  arise  in 
correlating  with  various  models  used  to  predict  a,  but  this  not  a  subject  to  concern  us 
here. 
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Figure  7. 
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LCF  data  for  60/40  solder  ( From  Ref.  2).  Plastic  strain  vs.  fatigue  life  (defined 
by  a  50%  drop  in  load). 


The  A 7p  -  Nf  curve  (the  correct  solid  curve  of  Figure  7)  can  be  converted  to  the 
required  A77  -  Nf  curve  by  adding  the  approximate  elastic  strains,  to  the  plastic 
strains,  A 77.  Care  must  be  taken  in  making  this  addition.  The  elastic  strain  is  depen¬ 
dent  upon  the  stress  level  which  is  in  turn  determined  by  the  deformation  rate  as  well 
as  the  amount  of  deformation.  In  fact,  solders  are  viscous  [4],  their  flow  stress 
depends  upon  the  strain  rate  and  not  just  the  strain  applied.  Fortunately  strain  rate 
data  was  generated  [4]  along  with  the  LCF  data.  Table  2  lists  the  shear  stress  range 

/\f 

Ar  required  to  produce  A 7*  =  — ,  where  G  was  calculated  from 

G 

^itrbr1'2*106  i6' 

with  v  =  0.5  (we  are  in  the  plastic  range)  and  E  =  4.5  x  10"6  -  1.5  x  106  psi  depending 
upon  the  reference  [16,17].  A  value  of  3.6  x  10"6  agrees  with  the  value  used  in  the 
FEA  and  it  yields  G  =  1.2  x  106  psi  which  is  close  to  the  1.6  x  106  psi  value  which 
correlated  the  data  of  Reference  2,  The  data  of  Table  3  was  used  to  construct  the  A77 
-  N f  curve  shown  in  Figure  8.  The  A77  curve  was  constructed  by  adding  the  four  N~je 
values  to  the  corresponding  A7 />  values  and  shifting  the  A77  curve  (the  a  =  0.59  curve 
was  used  because  it  used  the  correct  curve  fitting  procedure)  to  the  appropriate  A77 
level. 

Table  4  lists  Nf  for  the  strains  A77,  A772  and  A773  of  Table  1  (bottom  row). 

These  fatigue  lives  are  not  the  number  of  cycles  to  grow  a  fatigue  crack  through 
regions  1,  2  or  3.  This  is  because  there  is  an  important  size  effect  that  must  be  taken 
into  consideration. 

The  fatigue  life  of  Figures  7  and  8  was  defined  as  the  number  of  cycles  required  to 
reduce  by  1/2  the  load  required  to  produce  a  constant  strain  per  cycle,  A77.  The  load 
is  reduced  because  fatigue  cycling  results  in  the  nucleation  and  growth  of  fatigue 
cracks  which  reduce  the  load  bearing  area.  The  drop  in  load  represents  the  amount  of 
cracking.  A  criteria  of,  say,  a  25%  drop  in  load  would  correspond  to  less  cracking  than 
that  required  to  produce  a  50%  drop  in  load  and  fewer  fatigue  cycles  would  be 
required  for  only  a  25%  drop  in  load.  Likewise,  a  90%  drop  in  load  (i.e.,  when  the 
load  required  to  produce  A7 <p  is  only  1/10  of  that  initially  required)  would  necessitate 
more  cycles  than  that  for  a  50%  drop  in  load.  This  presents  a  problem  when  utilizing 
such  LCF  data.  What  load  drop,  or  other  parameter,  should  be  used  to  define  Nf  and 
what  definition  allows  for  the  best  application  of  this  data?  This  problem  is  not  a 
significant  one  when  other  materials  are  tested.  This  is  because  most  tests  are  run  in 
tension  on  materials  which  are  not  as  ductile  as  solder,  where  stable  (no  change  in  hys¬ 
teresis  load)  hysteresis  behavior  is  noted  over  a  large  number  of  cycles.  When  the 
load  begins  to  drop,  it  generally  does  so  relatively  rapidly  so  that  the  exact  definition 
used  to  define  Nf  does  not  appreciably  influence  the  data.  This  is  not  the  case  for 
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Figure  8.  Plastic  and  total  strains  vs.  Fatigue  life  (defined  by  a  50%  drop  in  load). 


Table  3 

Shear  Stress,  AE,  Tor  Various  Plastic  Shear  Strains,  Anp 
measured  at  35°C  and  03  Hz  (See  Ref.  2). 


Anp 

Ar  psi 

A nE 

A nr 

0.1 

10,000 

0.0083 

0.1083 

0.05 

9,000 

0.0075 

0.075 

0.02 

7,800 

0.0065 

0.0265 

0.01 

6,500 

0.0054 

0.0154 

Table  4 

ves  For  Strains  A~iru 

Anri  and 

Ann 

=  0.182 

Nf 

=  54 

AlT2 

=  0.095 

Nf- 

=  180 

Ann 

=  0.028 

II 

1900 
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solder  tested  in  shear  [2\.  In  the  shear  tests  performed  on  solder,  the  load  begins  to 
drop  after  the  first,  or  at  most  the  first  few,  cycles  and  continues  to  decline,  relatively 
slowly,  in  a  more  or  less  steady  manner.  The  number  of  cycles  required  to  drop  the 
load  by  25%  is  far  less  than  that  required  to  drop  the  load  by  90%. 

It  is  necessary  to  correlate  the  load  drop  definition  of  Nf  with  the  amount  of  crack¬ 
ing  required  to  produce  this  load  drop,  so  that  the  cycle  life  can  be  normalized  to 
account  for  the  different  amounts  of  cracking.  In  doing  this,  the  over-riding  assump¬ 
tion  is  that  all  the  load  drop  is  due  to  cracking  and  the  reduction  in  area  this  produces. 
We  are  thus  assuming  that  there  is  no  cyclic  hardening  or  softening  and  that  the  flow 
stress  remains  constant.  A  drop  in  load  can  thus  be  directly  related  to  a  decrease  in 
load  bearing  area  The  assumption  that  the  load  drop  was  due  to  crack  growth  was 
borne  out  in  a  qualitative  way  by  ultrasonic  microscopy  [2].  Figure  9  (taken  from 
Reference  2)  shows  an  example  of  the  results  of  an  ultrasonic  microscopy  of  a  sample 
which  was  cycled,  removed  for  examination,  recycled,  re-examined,  etc.  until  final 

AP 

failure.  The  parameter  <t>  =  1  -  — - - defines  the  load  drop.  It  is  equal  to  0  (i.e.,  0 

A/  max 

load  drop)  at  the  start  of  the  test  where  AP  =  APmax  and  1  (100%  load  drop)  when 
AP  =  O.  The  inserts  show  the  microscopy  results.  The  0.1  in.  x  0.5  in.  joint  area  is 
outlined  (with  the  shearing  in  the  0.5  in.  direction).  The  dark  areas  are  the  regions 
which  the  microscopy  delineates  as  cracked.  As  can  be  seen,  as  the  cycling  progresses, 
the  cracks  grow.  Unfortunately,  the  quality  and  resolution  of  the  ultrasonic  micro¬ 
graphs  does  not  permit  a  quantitative  analysis  (all  of  the  cracking  is  probably  not  being 
resolved). 


In  order  to  correlate  the  load  drop  <f>  with  a  crack  length,  it  is  necessary  to  make 
some  assumptions  and  simplifications  regarding  the  geometry  of  the  cracks.  If  it  is 
assumed  that  the  cracks  go  across  the  0.1  in.  joint  width  and  propagate  in  the  0.5  in. 
length  direction,  then  the  load  drop  can  be  correlated  with  a  crack  length  [2].  A  drop 
in  load  of  50%  ( <j>  =  0.5)  requires  a  crack  length  of  0.25  in.  (one-half  the  0.5  in  joint 
length)  if  a  single  crack  propagates  from  one  end  of  the  joint  or  0.125  in.  if  two  cracks 
are  present,  one  growing  from  each  end.  Figure  9  shows  that  both  of  these  assump¬ 
tions  are  simplifications  of  the  actual  behavior.  The  single  crack  assumption  is  better 
at  large  load  drops  where  the  cracks  have  coalesced  and  a  single  crack  may  dominate, 
but  at  very  small  load  drops  more  than  two  cracks  may  be  present. 


Nf  can  be  corrected  for  crack  growth  by  considering  the  average  crack  growth 
It  has  been  shown  [2]  that  as  inferred  from  the  drop  in  load  with  cycling, 

is  primarily  a  function  of  the  applied  strain  A%.  For  a  constant  A r/p,  is  thus  con 


AN 

stant.  (This  is  strictly  true  only  for  Ayp  at  25%  or  greater,  at  lower  strains 


A a 
AN 
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decreases  with  cycling.  This  effect  is  corrected  for  and  will  be  discussed  later).  For  O 
to  Nf  cycles,  i.e.,  AN  =  Nf,  and  a  crack  growth  of  D,  i.e.,  Aa  =  D. 


A a  _  D 
AN  ~  Nf 


[7] 


If  is  just  a  function  of  A~ip,  then  for  any  constant  A ~ip,  D  /Nf  can  be  related  to  the 
number  of  cycles,  Nj,  to  grow  a  crack  through  a  distance  L,  by 


Equation  8  can  be  rewritten  as 


D_  _  U. 
Nf  =  N, 


[8] 

[9] 


which  enables  one  to  calculate  Nj  from  L,,  D  and  Nf.  Equation  9  correlates  Ay  as 
measured  from  Figure  8  with  the  amount  of  cracking,  D,  used  to  define  Nf  and  the 
size  of  the  crack,  L,,  which  defines  Nj. 


Now  we  can  estimate  the  fatigue  life  for  the  process  depicted  in  Figure  6.  Let  us 
first  consider  the  cracking  under  the  Chip  Carrier  (i.e.,  region  1  in  Figure  6).  Mea¬ 
surements  of  the  CC/PWB  geometry  show  that  the  length  of  the  solder  under  the 
Chip  Carrier,  L  t,  is  0.04  in.  and  that  the  strain  in  this  region  A^r  =  0-182  in.  (see 
Table  2).  If  this  strain  were  to  have  been  applied  to  the  simple  solder  pads,  whose 
fatigue  behavior  is  described  in  Figures  7  and  8,  then  the  number  of  cycles,  Nf,  to 
reduce  the  load  by  50%  is  given  in  Table  4,  i.e.,  Nf  =  54  cycles.  The  number  of  cycles 
to  grow  the  crack  under  the  Chip  Carrier,  N  j,  from  Equation  9,  i.e., 


N  i 


0.04(54)  =  2A6 
D  D  ' 


Since  Nf  was  defined  by  a  50%  drop  in  load  for  an  0.5  in. 


long  solder  joint,  the  largest  value  for  D  is  0.250  in.  This  corresponds  to  a  single  crack 
growing  through  the  pad.  Figure  9  shows  that  D  should  be  smaller  than  0.250  in, 

reflecting  the  existence  of  multiple  cracks.  For  four  cracks,  D  =  -  —  =  0.0625 


in.,  which  represents  a  D  value  which  should  be  at  the  low  end  of  the  probable  values 
for  D  for  cracks  which  produce  a  50%  drop  in  load.  Multiple  cracking  is  even  more 
pronounced  at  the  start  of  the  test,  but  less  so  at  the  end  of  the  test  where  more  crack 
coalescence  has  taken  place.  The  number  of  cycles  to  grow  the  crack  into  the  fillet, 

N 2,  and  through  the  fillet  N 3  can  be  calculated  in  a  similar  manner.  These  results, 
with  D= 0.062  in.  and  D= 0.250  in.,  are  summarized  in  Table  5.  The  lengths  Li  and 
L3  were  estimated  from  the  Chip  Carrier  geometry  and  the  FEA  strain  distribution 
(which  was  used  to  separate  the  transition  region  from  the  rest  of  the  fillet).  The  Nf 
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Correlation  of  ultrasonically  measured  cracked  areas  and  the  drop  in  load 
measured  in  interrupted  tests. 


Table  5 

Calculation  of  A',  From  Nr,  L,  and  L) 


values  were  determined  from  the  A 77,-  values  as  shown  in  Table  4.  A/,  not  Nf  is  corre¬ 
lated  with  the  results  of  fatigue  cycling  actual  joints. 


As  has  already  been  noted,  — —  is  not  only  a  function  of  A-yp.  At  strains  below 

AN 

A/i 

25%,  ,  as  determined  from  the  load  drop,  decreases  with  increasing  number  of 

cycles.  This  gives  rise  to  differences  in  the  A-t />  -  Nf  curves  when  different  choices  are 
made  for  the  definition  of  Nf.  Figure  10  shows  such  curves  for  Nf  defined  by  a  25%, 
50%  and  90%  drop  in  load  criteria  (the  data  for  the  50%  curve  is  shown  in  Figure  7). 
As  can  be  seen,  a  decreases  when  the  load  drop  definition  for  Nf  increases.  This 
occurs  because  the  lower  the  strain,  A-yp,  the  more  the  rate  of  load  drop  decreases 

with  cycling,  producing  a  decrease  in  (A  separate  publication  is  currently  being 

prepared  to  describe  the  relationship  of  a  and  the  degree  of  load  drop  used  to  define 
Nf).  The  Ayr  curves,  created  from  the  calculated  Aye  are  also  shown  in  Figure  10  and 
can  be  used  to  calculate  Nf  for  a  25%,  50%  and  90%  drop  in  load.  The  results  are 
listed  in  Table  6  for  the  same  A77  values  used  in  Table  4.  N  UN2  and  /V3  can  likewise 
be  calculated  in  the  same  way  as  before,  i.e.,  by  Equation  9.  The  results  of  such  calcu¬ 
lations  are  listed  in  Table  7.  There  is  good  agreement  between  the  Nl,  N2  and  A/ 3 
values  derived  from  the  25%,  50%  and  90%  Nf  -  A-yp  curves  of  Figure  10,  providing 
the  correct  value  is  assigned  to  D. 


A  single  crack  across  the  joint  width  would  have  to  grow  0.125  in  to  produce  a  25% 
drop  in  load.  As  noted,  however,  in  the  initial  stages  there  are  many  cracks  growing  so 
0.125/5  or  0.021  in  was  used  to  define  D.  By  the  time  the  crack  has  grown  to  reduce 
the  load  by  50%,  there  has  been  a  coalescence  of  the  cracks.  It  is  thus  assumed  that 

025 

three  cracks  have  grown  through  the  0.25  half-length,  so  D  -  — y—  =  0.083  in.  There 

is  further  coalescence  as  the  cracks  grow,  so  that  by  the  time  a  90%  drop  in  load  is 
reached,  there  is  one  dominant  crack.  It  is  thus  assumed  that  1-1/2  cracks  are  grow- 

0  45 

ing  through  a  0.45  in.  length,  so  D  =  -yyy  =  0.30  in.  The  choice  of  an  increasing  D 

value  with  load  drop  and  the  reduction  in  slope  of  the  Ay  -  Nf  curves  with  increasing 
load  drop  corrects  for  the  error  introduced  by  assuming  Aa/A N  was  only  a  function  of 
A-y p.  If  this  were  the  case,  the  A'y  -  Nf  curves  should  be  parallel  and  spaced  according 
to  the  load  drop  used  to  define  Nf  (i.e.,  the  50%  load  drop  curve  would  be  positioned 
at  twice  the  number  of  cycles  as  the  25%  curve).  The  N,  values  listed  in  Table  7  will 
be  compared  to  the  fatigue  lives  measured  in  tests  of  actual  CC/PWB  joints. 
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/^(L,  =0.04  m.)  Nz(L2=0m  in.)  N,(L}=  0.015  in 


For  Nf  (90%) 
D  =  0.30 


15 


15 


300 


III.  Comparison  with  CC/PWB  Joint  Tests 

The  main  reason  why  the  preceding  analysis  was  performed  was  to  predict  the 
behavior  of  CC/PWB  joint.  Such  joints  have  been  subjected  to  LCF  by  applying  shear 
displacements  across  the  joints.  The  results  of  these  tests  have  been  discussed  exten¬ 
sively  [13].  Furthermore,  the  joints  were  made  with  approximately  a  3-mil  height  of 
the  CC  above  the  PWB,  which  is  what  was  assumed  for  the  FEA  calculations.  The 
actual  fillet  geometry  was  between  the  bulbus  and  nominal  geometries,  hence  the 
averaging  of  the  calculated  strains  that  was  used  previously. 

There  was  a  difference  in  the  solder  and  microstructure  between  that  of  the  LCF 
data  that  was  used  in  the  calculation  and  what  was  used  in  the  actual  CC/PWB  joints. 
The  previous  LCF  data  utilized  60Sn/40Pb  solder  with  a  coarse  microstructure.  The 
CC/PWB  joint  tests  utilized  63Sn/37Pb  with  a  finer  microstructure.  Numerous  tests 
performed  in  the  rest  of  the  Man-Tech  program  have  shown,  however,  that  these  are 
not  significant  effects  insofar  as  the  35°C  isothermal  fatigue  behavior  is  concerned 
(this  is  what  we  are  trying  to  correlate).  Hence,  we  shall  use  the  60/40  LCF  data  for 
predictions,  without  any  correction  for  the  slight  difference  in  composition  and  micro¬ 
structure  between  this  and  that  of  the  CC/PWB  joints. 

In  this  section  we  shall  compare  the  predictions  of  joint  lives  made  previously  with 
the  lives  actually  measured  in  the  tests  on  CC/PWB  joints.  The  common  thread  is 
that  the  correlations  will  be  made  for  the  same  displacements.  Before  the  test  results 
are  considered,  it  is  necessary  to  correct  for  the  fact  that  the  CC/PWB  tests  were  run 
with  a  constant  plastic  displacement  A />,  whereas  the  FEA  calculations  were  done  for 
a  given  total  displacement  A  r- 

The  detailed  specifics  given  below  refer  to  the  specifics  of  the  experiments  that 
were  run  and  the  FEA  calculations  that  were  made.  They  may,  therefore,  not  be  gen¬ 
erally  applicable.  They  have  been  included,  however,  not  only  to  illustrate  exactly  how 
the  life  calculations  were  made,  and  any  approximations  that  were  made,  but  also  to 
emphasize  the  importance  of  making  sure  that  the  experimental  conditions  and  FEA 
assumptions  coincide  or  that  any  differences  are  corrected  for. 

The  conversion  of  A p  to  At  is  complicated  by  the  fact  that  when  Ap  is  kept  con¬ 
stant  throughout  the  testing,  as  was  approximately  done  here,  the  elastic  displacement 
decreases  due  to  the  drop  in  load  which  accompanies  the  cycling.  It  is  therefore 
necessary  to  determine  the  relationship  of  A p  and  A7  at  various  stages  during  the  test 
and  correlate  this  with  the  number  of  fatigue  cycles  reached  by  that  point  in  the  test. 

The  elastic  displacement  is  determined  by  the  elastic  compliance  K0,  i.e., 

Ae  =  K0(£Joad)  [10] 

This  compliance  contains  displacements  of  the  solder,  PWB,  CC  and  any  bending  that 
was  present.  Measurements  [13]  were  made  on  an  uncut  PWB  with  no  CC  or  solder 
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joint.  These  measurements  and  calculations  of  the  elastic  behavior  of  the  PWB  and 
CC  yielded  the  result  that  the  compliance  due  to  the  solder  alone  was  given  by 

tfs=0.68 K0  [11] 

so  the  elastic  compliance  due  to  the  solder  A Es  is 

A ES  =  0.68*o  (A/oad)  [12] 

Since  all  the  plastic  flow  is  due  to  the  solder  (the  tests  on  the  uncut  PWB  showed  that 
f  for  loads  in  excess  of  what  was  actually  used,  the  PWB  behaved  elastically)  the  total 

n  displacement  used  in  the  FEA  can  be  related  to  the  measured  plastic  displacement 

and  load  range  by: 

,  At  =  AP  +  0.68*o  (Aload)  [13] 

.  Figure  11  shows  the  load  range  A(load)  as  a  function  of  the  plastic  strain  range, 

^  both  of  which  were  measured  at  the  start  of  the  test.  The  data  was  taken  at  the  start 

of  each  experiment  and  represents  the  behavior  of  each  test  sample.  The  displace¬ 
ment  A p  was  measured  across  both  rows  of  joints.  This  is  the  displacement  developed 
by  thermal  mismatch  and  the  displacement  utilized  in  most  experiments.  It  is  not, 
however,  the  displacement  utilized  in  the  FEA.  There  the  displacement,  APS,  across  a 
single  row  of  joints,  was  considered.  At  the  start  of  the  test,  in  the  symmetrically 
'  prepared  joints  utilized  here,  the  applied  displacement  was  more  or  less  equally  dis- 

;  tributed  between  the  two  rows  of  joints  [  13].  By  using  an  arrangement  which  allowed 

the  measurement  of  the  displacements  across  each  row  separately  it  was  found  that 

A  PS  =  0.49Ap  [14] 

where  A p$  is  the  plastic  strain  measured  across  the  single  row  of  joints  which  eventu¬ 
ally  fails.  The  compliance  measured  across  a  single  row  of  joints,  at  the  start  of  a  test, 
i  is  also  about  one-half  that  measured  across  both  rows  of  joints.  The  elastic  displace¬ 

ment  of  the  solder  across  a  single  row,  AEss  is  thus  (using  a  0.49  correction  factor) 


A  Ess  =  0.33KoA(load)  [15] 

The  total  displacement,  Aps,  measured  across  a  single  row  of  joints  is  therefore 

A ts  ~  &rs  +  &ess  [16] 

A  ts  =  A  PS  +  0.33KoA(load)  [17] 

Ats  =  0.49Ap  +  0.33KoA(load )  [18] 


Equation  16  relates  A 75,  which  was  used  in  the  FEA  with  the  plastic  and  elastic 
displacements  of  the  solder  of  a  single  row  of  joints.  Equations  17  and  18  relates  A  73 
with  the  compliance  measured  across  both  rows  of  joints,  the  plastic  strain  measured 
across  both  rows  of  joints  and  the  load  range  of  the  initial  hysteresis  loop.  The  plastic 
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Figure  11.  Load  range  of  the  initial  hysteresis  loops  vs.  applied  plastic  displacement  AP. 
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strain  Ap  and  the  ioad  range  A(load)  are  correlated  by  the  data  of  Figure  1 1.  The 
best  fit  to  this  data  yields: 

A  (load)  =  757.4A22007  [19] 

or  for  the  measured  [13]  value  of  K0  ~  7.20  x  10-6  in/lb. 

Ary  =  0.49Ap  +  1.80  x  10 ^A^2007  [20] 

Equation  20  relates  A ts  (used  in  the  FEA)  with  A p  (used  as  the  control  limit  for  the 
tests  on  the  CC/PWB  joints). 

The  constants  were  derived  for  the  start  of  the  test.  They  change,  however,  with 
continued  cycling.  As  the  cycling  continues,  the  joints  crack.  This  causes  the  load 
required  to  produce  a  given  Ap  to  decrease  and  the  compliance  to  increase.  Because 
much  of  the  measured  compliance  is  not  due  to  the  solder  deflections  and  generally 
only  row  of  joints  cracks  (resulting  in  a  compliance  change  only  for  that  row),  the 
compliance  change  is  much  less  than  would  be  expected  from  the  load  drop  (i.e.,  when 
the  load  drops  in  half  the  compliance  does  not  double). 

The  decrease  in  compliance  causes  the  applied  Ap  to  decrease  slightly.  The  tests 
were  run  with  Ap  limits  with  Ap  being  determined  by  the  initial  compliance.  The 
change  in  compliance  results  in  not  enough  elastic  strain  being  subtracted  from  the 
total  strain  signal  and  a  decrease  in  the  control  Ap  (see  Reference  13  for  more  discus¬ 
sion  of  this  effect).  Figure  12  shows  the  increase  in  K  (as  measured  across  both  joints) 
and  the  Ap  decrease  resulting  from  this  K  increase.  Also  shown  is  Ks  the  compliance 
measured  across  the  joint  that  failed.  Because  the  cracking  is  concentrated  in  this 
row,  Ks  shows  a  larger  increase  than  K.  Also  shown  is  Apy/AP  which  shows  how  the 
plastic  strain  becomes  concentrated  in  the  row  of  joints  which  crack.  All  the  variables 
(except  Aps/Ap  )  are  shown  relative  to  the  initial  values  measured  at  the  start  of  the 
test.  Figure  12  shows  how  K /Ka,  Ks /Kso,  Ap/Apo  and  Apy /Ap  vary  with  cycling  as 
defined  by  the  drop  in  load.  There  is,  however,  relatively  little  change  in  any  of  these 
variables  below  about  a  50%  drop  in  load. 


The  data  of  Figure  12  can  be  utilized  to  correlate  A 7-$  with  Ap  or  Apy  for  any  load 
drop.  Equation  20  can  be  written  as  a  generalized  expression 


A  ts  - 


Apy' 

Ap  + 

Apy 

0.68 

K) 

,  , 

K0\{load) 


[21] 


where  the  expressions  in  the  brackets  are  described  in  Figure  12.  One  further  correc¬ 
tion  is  required  because  A(Ioad)  is  in  general  not  given  entirely  by  Equation  19.  This 
expression  describes  the  load  at  the  beginning  of  the  test.  For  any  load  drop  LD,i t 

must  be  corrected  by  a  ( factor,  (i.e.,  the  load  is  reduced  by  0.75  for  a  25%  load 

lw 

drop).  Substituting  this  factor  and  Equation  19  into  Equation  21,  one  gets 
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Figure  12.  Changes  in  the  compliance,  compliance  across  the  joint  that  is  cracking,  plastic 
strain  range  and  plastic  strain  across  the  joint  that  is  failing  vs.  the  %  drop  in 
load. 
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[22] 


Equation  22  is  a  generalized  expression  which  relates  A ps  with  A />  for  any  load 
drop  (Equation  20  is  valid  only  for  the  start  of  the  test  where  Aps/Ap  -  0.49,  —  = 


1,  and  LD  -  O).  This  equation  enables  us  to  relate  the  FEA  performed  for  a  given 
A ts  with  the  CC/PWB  tests  performed  with  A p  limits  for  any  point  in  the  test  (i.e., 
for  any  drop  in  load).  The  fatigue  life  was  measured  for  25%,  50%  and  90%  drop  in 
load. 


Figures  13  and  14  show  the  results  for  Nf  defined  by  a  25%  drop  in  load.  For  this 
drop  in  load,  K /K0  =  Ks/Kso  =  1,  Ap/Apo  =  l,  and  Aps/Ap~ 0.49.  Only  the  load 
drop  term  is  appreciably  different  from  its  initial  value.  It  decreases  to  0.75  from  1,  so 
an  average  value  of  0.875  will  be  used  in  Equation  22,  i.e.,  with  K0  =  7.20  x  104*  in/lb, 
Equation  22  becomes 

A  TS  =  0.49Ap  +  1.59  x  10'3a£2007  [23] 


For  Ajs  =  5  x  10"4  in.  (twice  the  value  used  in  the  FEA  which  considered  only  0  to  + 
A  whereas  the  experiments  were  performed  for  -A  to  +  A  ).  A p  =  3.60  x  10-4  in.  and 
ApS  =  1.76  x  1C4  in.  These  are  the  plastic  displacements  which  correspond  to  an 
applied  total  displacement  across  one  row  of  joints  of  5  x  10"4  in. 

Figures  13  and  14  can  be  used  to  determine  Nf  for  these  plastic  displacements. 

For  Nf  defined  by  a  25%  drop  in  load,  these  curves  yield  Nf  =  26  for  the  AP  data  (i.e., 
with  the  plastic  displacement  measured  across  both  rows  of  joints)  and  Nf  =  33  cycles 
for  the  correlation  with  the  displacements  measured  across  the  row  of  joints  which 
actually  failed.  The  differences  between  these  two  measurements  reflects  differences 
in  the  curve  fits  to  the  data.  The  average  value  of  Aps /A p  =  0.49  but  the  data  curves 
are  displaced  by  a  factor  of  0.53.  The  Aps  correlation  is  most  appropriate  for  compar¬ 
isons  with  the  FEA  which  was  done  on  a  single  joint.  The  correlation  with  A p  is  given 
because  this  is  what  would  be  calculated  in  a  thermal  displacement  calculation.  Ini¬ 
tially  the  displacement  is  equally  shared  by  both  rows  of  joints  so  Ap~2APS.  As  the 
test  progresses,  however,  more  of  the  displacement  becomes  shared  by  the  row  of 
joints  in  which  cracks  are  developing. 

An  examination  [13]  of  specimens  has  shown  that  a  load  drop  of  25%  correlates 
with  cracking  under  the  chip  carrier,  i.e.,  with  region  1  of  Figure  6.  These  Nf  values 
should  therefore  correlate  with  N  {  from  Table  7  (Column  1).  The  measured  value  of 
26-33  cycles  is  in  excellent  agreement  with  the  predicted  values  of  15-38  cycles  (see 
Table  7). 
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•e  13.  Plastic  displacement,  AP,  across  both  rows  of  joints  vs.  Nf  for  a  25%  drop  in 
load. 


Figure  14.  Plastic  displacement,  APs,  across  the  row  of  joints  that  failed  vs.  Nf  for  a  25% 
drop  in  load. 
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A  drop  in  load  of  50%  resulted  in  cracks  which  grew  into  the  CC/PWB  joint  fillet 
but  did  not  propagate  through  the  fillet.  We  will  therefore  correlate  this  behavior  with 
the  number  of  cycles  to  grow  through  regions  1  and  2,  i.e.,  Nf  for  a  50%  drop  in  load 
should  correlate  with  N i  +  A2  as  listed  in  Table  7.  When  the  load  drops  by  50%, 
there  are  significant  changes  mK/K0,  Ks/Kso,  Ap/Apo  and  Apj/Ap  as  listed  in  Fig¬ 
ure  12.  The  average  changes  from  0  to  a  50%  drop  in  load  are  K /K0  -  1.02,  Ks  /Kso 

=  1.06,  Ap/Apo  =  0.924,  and  ApS/AP  =  0.525,  [  1-^ 


100 


=  0.75  (the  average  LD  is 


25%).  Equation  22  thus  reduces  to 


A ts  =  0.525A p  +  1.489  x  lO^2007  [24] 

For  Ars  =  5  x  10"4,  A p  =  3.71  x  10-4  and  with  A P$  =  0.525A p  =  1.95  x  10"4.  Fig¬ 
ures  15  and  16  plot  Nf  for  a  50%  drop  in  load  in  the  CC/PWB  joint  tests  vs.  A p  and 
A ps-  Here  the  plastic  displacements  are  the  averages  for  zero  drop  in  load  to  50% 
drop  in  load.  For  A p  =  3.71  x  10"4  in.,  Nf  =  64  cycles  and  for  A.PS  =  1.95  x  10-4  in., 

Nf  =  59  cycles.  These  results  are  also  in  excellent  agreement  with  N{  +  N2  from 
Table  7,  i.e.,  with  30-64  cycles. 


By  the  time  the  load  has  dropped  by  90%,  cracks  are  visible  in  most  of  the  fillets  in 
the  CC/PWB  joints  [14].  This  load  drop  should  thus  correlate  with  A7  j  +  N2  +  N3  of 
Table  7.  The  rate  in  the  drop  in  load  slows  down  after  a  50%  drop  in  load  is  reached 
[14],  It  takes  many  more  cycles  to  reduce  the  load  from  50%  to  90%  than  to  reduce  it 
from  0  to  50%.  Thus,  instead  of  considering  the  average  value  of  the  parameters  of 
Figure  12  over  0  to  a  90%  drop  in  load,  we  will  use  the  average  from  50%  to  90% 
drop  in  load  and  likewise  Figures  17  and  18  show  Nf  defined  by  a  90%  drop  in  load  as 

i  .  ju  ^p(.5)  +  Ap(.9)  Aps(.5)  +  Ap$(.9) 

correlated  by - - - and - - - .  For  K /K0  =  1.305,  KsJKso 

=  1.369,  Ap/Ap0  =  0.806,  APS/AP  =  0.655  and 
becomes 


LD 

100 


=  0.3,  Equation  22 


A  ps  =  0.655Ap  +  9.509  x  lO^A^2007  [25] 

For  Ars  =  5  x  10"4,  AP  =  4.54  x  104  and  for  APS  =  0.655,  AP  =  2.97  x  104.  Because 
of  the  drop  in  load  and  concentration  of  the  displacement  in  the  row  of  joints  that  is 
failing,  A Ps  is  approaching  Ap  and  AP  is  approaching  Aps-  For  AP  =  4.54  x  10'4  in., 

Nf  =  285  cycles  and  A p$  =  2.97  x  1CT4  in.,  Nf  =  280  cycles.  Because  of  hysteresis 
loop  distortions  which  developed  at  large  load  drops,  Nf  was  also  correlated  with  Ap 
as  measured  at  zero  load  instead  of  at  the  load  reversal  (see  Reference  13).  The  Nf 
(90%  drop  in  load)  as  correlated  by  this  zero  load,  Ap  measurement  is  shown  in  Fig¬ 
ure  19.  For  Ap  =  4.54  x  104  in.,  Nf  -  200  cycles.  These  Nf  values  are  in  good  agree¬ 
ment  with  N  ]  +  N2  +  Af  3  as  listed  in  Table  7  (i.e.,  313  to  390  cycles). 


26 


200 

100 


* 


} 


* 


Figure  15.  Plastic  displacement.  A,  across  both  rows  of  joints  vs.  Nf  for  a  50%  drop  in 
load. 


Figure  16.  Plastic  displacement,  Aps,  across  the  row  of  joints  that  failed  vs.  Nf  for  a  50% 
drop  in  load. 
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Figure  17.  Plastic  displacement.  A,  across  both  rows  of  joints  vs.  Nr  for  a  90%  drop  in 
load. 
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Figure  18.  Plastic  displacement,  Are  across  the  row  of  joints  that  failed  vs.  Nf  for  a  90% 
drop  in  load. 
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Figure  20  compares  all  the  predictions  and  measurements,  and  illustrates  the 
excellent  agreement  that  has  been  obtained  with  predictions  falling  within  ±2X  of  the 
measured  values.  The  measured  data  was  obtained  from  the  A p  and  A/>s  correlations 
and  the  small  differences  represent  curve  fit  differences  to  the  data.  It  is  important  to 


remember  that  A p$  is  not  just 


Ap  A P 

— -  and  if  — — 
2  2 


were  used  to  define  the  behavior  of  a 


single  row  of  joints  (the  condition  for  which  the  calculations  were  made),  a  serious 
error  would  be  made.  If  only  A/>  data  were  available,  some  sort  of  estimate  would 
have  to  be  made  to  determine  how  the  strains  were  partitioned.  Equal  partitioning 
may  be  o.k.  at  the  start  of  the  test,  but  it  is  not  a  correct  procedure  to  follow  at  later 
lives.  The  most  conservative  approach  would  be  to  assume  A p  =  A/>s  which  would 
yield  an  underestimate  of  the  actual  fatigue  life. 


The  values  of  D  used  in  Table  7  greatly  influence  the  predicted  fatigue  lives.  The 
D  values  were  chosen  to  give  the  best  balance  of  predicted  lives  vs.  the  measured  ones. 
This  parameter  is  not,  however,  just  an  error  minimization  factor.  It  must  conform  to 
the  size  of  the  pad  used  in  generating  the  Nf  data  and  to  the  degree  of  multiple  crack¬ 
ing  at  any  stage  in  life.  The  accuracy  of  the  predictions  also  depends  upon  the  accu¬ 
racy  with  which  the  FEA  determines  the  strain  distribution  and  the  degree  to  which 
the  joint  geometry  can  be  specified  accurately.  Even  if  a  very  accurate  prediction  can 
be  made,  at  best  it  will  only  reflect  the  mean  expected  behavior.  The  scatter  in  the 
experimental  data  is  such  that  the  behavior  of  individual  specimens  can  significantly 
deviate  from  this  mean. 


The  scatter  in  the  measured  values  of  Nf  can  be  greater  than  3X,  depending  upon 
the  confidence  limits  that  are  considered.  Figures  13-18  show  the  least  squares  best  fit 
trend  line  and  dotted  lines  which  define  70%  and  95%  confidence  limits.  These  are 
close  to  but  not  exactly  the  same  as  one  and  two  standard  deviation  limits.  A  70% 
confidence  boundary  means  that  there  is  70%  confidence  that  the  results  for  the  next 
test  will  fall  within  this  band  and  a  15%  probability  that  it  would  fall  above  the  bound¬ 
ary  and  a  15%  probability  that  it  would  fall  below  the  boundary.  Likewise,  for  the 
95%  confidence  zone,  there  is  a  2.5%  probability  that  the  next  data  point  would  be 
below  this  band.  In  the  data  sets  with  the  most  scatter,  there  is  a  factor  of  5  spread 
between  this  95%  confidence  boundary  and  the  mean.  Thus,  there  is  a  2.5%  probabil¬ 
ity  that  a  measured  value  could  be  only  1/5  the  mean  value.  This  is  not  a  particularly 
surprising  fact,  nor  is  the  level  of  scatter  particularly  large  considering  the  real  joints 
are  being  tested,  not  carefully  machined  specimens.  A  design  which  calls  for  greater 
than  97.5%  confidence  requires  more  than  a  factor  of  safety  of  5X  on  measured  mean 
values. 
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Figure  19.  Plastic  displacement,  Ap2,  across  both  rows  of  joints  as  measured  at  zero  load 

vs.  Nf  for  a  90%  drop  in  load. 


Figure  20.  Comparison  between  the  predictions  made  from  the  FEA  and  LCF  data  and 
measurements  on  CC/PYVB  joints. 
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The  scatter  in  the  data  makes  it  difficult  to  detect  subtle  effects  such  as  slight 
changes  in  solder  composition.  Gross  changes  such  as  large  changes  in  the  fillet 
geometry  can,  however,  be  measured.  These  effects  will  be  discussed  in  subsequent 
Man-Tech  papers. 


IV.  Predicting  Thermal  Fatigue  Lives 

All  the  aforementioned  data  and  discussion  refer  to  tests  run  at  one  temperature 
where  the  fatigue  results  from  mechanical  displacements.  While  a  knowledge  of  this 
behavior  is  critical  to  predicting  the  thermal  fatigue  behavior  of  CC/PWB  joints,  it  is 
not  enough  to  make  accurate  predictions.  In  this  section  we  shall  discuss  the  method¬ 
ology  for  making  predictions  of  the  thermal  fatigue  of  joints  and  some  simple  approxi¬ 
mations  that,  together  with  the  isothermal  approaches  described  in  this  paper,  could 
be  used  to  make  a  first,  rough,  estimate  of  the  fatigue  life.  The  steps  required  to 
make  a  thermal  fatigue  life  estimate  are  as  follows: 


A  Determine  the  Thermal  Displacement 


The  previous  analysis  has  shown  that  in  order  to  estimate  the  fatigue  life  it  is 
necessary  to  know  the  displacement  being  applied  to  a  solder  joint.  This  is  as 
true  for  thermal  fatigue  as  it  is  for  the  isothermal  fatigue  considered  here.  This 
thermal  displacement  arises  because  of  the  mismatch  in  the  coefficients  of  ther¬ 
mal  expansion,  Aa,  between  the  CC  and  PWB.  In  a  leadless  device,  this  mis¬ 
match  is  taken  up  by  the  solder  joint  and  deflections  in  the  PWB.  There  is  no 
lead  to  deflect  and  take  up  this  displacement.  To  a  first  approximation  for  a 
temperature  change  AT  the  thermal  displacement  can  be  calculated  as 


A  = 


(Aa)  (A T)L0 
P 


[26] 


where  La  is  a  CC  length  dimension,  typically  the  longest  distance  between 
solder  joints  (close  to  the  diagonal  length),  and  P  is  partition  coefficient  equal 
to  2  if  the  displacement  is  equal  on  opposite  sides  of  the  chip  carrier  and  equal 
to  1  if  all  the  displacement  is  assumed  to  be  on  one  side.  The  largest  plastic  dis¬ 
placement  A p  is  achieved  if  all  the  displacement  is  plastic  solder  joint  displace¬ 
ment  and  using  P  =  1 

A P  =  (Aa)  (AT)  L0  [27] 

Equation  27  assumes  that  all  the  displacement's  taken  up  on  one  row  of 
joints  and  that  all  of  these  displacements  are  plastic.  Above  room  temperature, 
solder  creep  relieves  any  thermal  stresses  and  the  assumption  that  all  the  dis¬ 
placement  is  plastic  is  a  good  one.  Below  room  temperature,  however,  there  is 
much  less  creep  and,  depending  upon  the  exact  temperature  and  times  involved 
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in  the  thermal  cycling,  there  may  be  significant  elastic  solder  strains  and 
deflections  of  the  PWB.  As  an  approximation,  it  can  be  assumed  that  all  the 
fraction  of  A  T  spent  above  room  temperature,  Tr,  produces  plastic  strains  but 
that  only  1/5  of  the  thermal  displacement  developed  below  Tr  produces  plastic 
displacements,  i.e., 


A/.  = 


Tr  -  7#j  AaATL0 

B  AT  5 

/ 


[28] 


(where  Tm  is  the  maximum  temperature  and  TN  is  the  minimum  temperature). 
For  equal  excursions  above  and  below  Tr,  i.e.,  Tm  -  Tr  =  Tr  -  Tn  =  ~~ 


Ap  = 


1 
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AaATL0 

5 


[29] 


Furthermore,  if  B  =  2  (i.e.,  half  the  displacement  developed  below  Tr  produces 
plastic  displacement)  and  P  =  1.5,  which  assumes  some  unequal  displacement 
partition  (between  5  =  1  which  denoted  partition  to  one  side  only  and  P  -  2 
which  is  equal  partition)  Equation  29  becomes 


A P  =  0.5(Aa)(AT)Lo  [30] 

which  is  a  reasonable  estimate  for  A p  (for  equal  thermal  cycles  above  and  below 
Tr),  but  one  which  may  still  overestimate  A p.  The  most  conservative  (largest) 
simple  estimate  for  AP  is  given  by  P  =  1  and  5  =  1  and  is  given  by  Equation  27 
no  matter  how  AT  is  distributed  above  or  below  Tr. 

As  a  first  estimate  of  AP  one  should  use  Equation  26  or  28.  If  the  predicted 
life  plus  a  safety  factor  exceeds  the  design  guidelines,  then  no  further  calculation 
of  A p  is  needed.  If  this  is  not  the  case,  then  a  more  refined  calculation  is 
needed.  Equation  29  may  overestimate  AP  because  P  may  be  closer  to  2  and  5 
closer  to  oo.  For  5  =  2  and  5  =  oo.  Equation  29  becomes 


A P  =  0.25(Aa)(Ar)Eo  [31] 

which  is  a  simple  lower  bound  estimate  for  A p  but  not  the  lowest  possible  value 
for  A p  because  it  neglects  the  possibility  that  some  of  the  thermal  strains  are 
taken  up  by  bending  the  PWB. 

Determining  the  exact  value  for  A P,  as  opposed  to  the  simple  approxima¬ 
tions  discussed  here  is  difficult,  but  can  be  done  [8,9,12].  A  proper  estimate  of 
A p  is  critical  because  Nf  varies  approximately  with  the  square  of  AP.  Thus, 
there  would  be  a  factor  of  16  difference  in  Nf  as  estimated  from  Equation  27  vs. 
Equation  31.  A  refined  calculation  of  a  whole  circuit  board  populated  by  a  large 
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number  of  leadless  chip  carriers  and  other  components  might  show  that  even 
Equation  27  was  not  an  upper  bound  to  A />  because  PWB  bending  might  be 
concentrated  in  a  region  giving  an  extra  contribution  to  A />  not  accounted  for  in 
Equation  27.  There  may  also  have  to  be  special  consideration  given  to  the 
comer  joints  as  these  may  also  experience  additional  displacement  due  to  asym¬ 
metries  in  bending,  heat  transfer,  and  the  largest  La  for  the  Chip  Carrier  (i.e., 
the  diagonal  length  rather  than  the  opposite  side  length). 

The  PWB  bending  causes  the  chip  carrier  to  experience  tensile  as  well  as 
shear  strains,  and  the  chip  tends  to  move  up  away  from  the  PWB.  This  compli¬ 
cated  displacement  pattern  should  serve  as  the  boundary  condition  for  the  calcu¬ 
lations  discussed  below. 

B.  Calculation  of  Strains  From  the  Thermal  Displacements 

Once  the  thermal  displacement  is  determined,  the  joint  strain  distribution 
can  be  determined  by  FEA.  This  calculation  is,  however,  considerably  more 
difficult  than  the  one  discussed  here.  All  the  elastic  and  plastic  properties  vary 
with  temperature.  Furthermore,  the  yielding  of  the  joint  will  be  a  function  of 
not  only  the  temperature  but  also  the  rate  of  temperature  change,  which  estab¬ 
lishes  a  strain  rate. 

This  calculation  should  be  due  to  yield  A  ip  rather  than  A77-.  In  the  calcula¬ 
tions  performed  here,  it  was  a  relatively  simple  manner  to  convert  A^P  -  Nf 
curves  to  An  -  Nf  curves.  When  the  temperature  varies,  so  does  the  elastic 
strain  (due  to  differences  in  the  yield  strength  and  elastic  modulus).  Thus  the 
conversion  of  A ~fP  data  to  A77  data  is  difficult  and  should  therefore  be  done  as 
part  of  the  FEA,  where  a  temperature  variation  can  be  incorporated. 

The  most  accurate  calculation  should  employ  the  complicated  displace¬ 
ments,  developed  because  of  PWB  bending,  as  the  boundary  conditions  for  the 
calculation.  Furthermore,  these  calculations  should  also  take  any  metallurgical 
instabilities  into  account.  In  thermal  fatigue  tests,  unlike  the  isothermal  tests 
discussed  here,  a  grain  coarsened  or  recrystallized  zone  often  develops  under 
the  Chip  Carrier  and  extends  into  the  fillet.  The  mechanical  properties  of  this 
zone  are  different  from  those  of  the  rest  of  the  joint.  Strains  become  concen¬ 
trated  in  this  zone  and  are  higher  than  would  be  estimated  were  no  local  metal- 
lurgically  distinct  zone  formed.  These  increased  strains  will  reduce  the  fatigue 
life.  Thus  the  formation  of  a  grain  coarsened  or  recrystallized  zone  is  an  impor¬ 
tant  factor  in  determining  the  fatigue  life,  and  is  a  factor  which  an  accurate  life 
prediction  must  account  for. 

If  the  mechanical  properties  of  the  coarsened  or  recrystallized  zone  can  be 
determined,  then  the  FEA  can  be  modified  to  account  for  this  zone.  The  FEA 
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can  then  yield  the  strain  in  this  one  and  these  strains  can  be  used,  following  pro¬ 
cedure  C  below,  to  calculate  Nf. 

C.  Calculation  of  Nf  from  A^P 

Once  A'lp  is  calculated  for  the  thermal  cycle,  it  is  a  relatively  simple  matter 
to  calculate  Nf,  following  the  procedures  laid  out  in  this  report.  A  distribution 
of  A^p  strains  should  be  considered  and  N;  calculated  for  each  relevant  strain  in 
the  distribution.  Fortunately  there  does  not  appear  to  be  a  large  influence  of 
temperature  on  the  A^p  -  Nf  behavior  [5],  and  to  a  first  approximation,  room 
temperature  data  can  describe  the  behavior  over  a  range  of  temperatures. 

There  is,  however,  a  correction  which  must  be  made  to  account  for  the  cycling  fre¬ 
quency.  In  a  typical  thermal  cycle  test  this  is  determined  by  the  rate  of  tempera¬ 
ture  change  and  hold  time  at  temperature.  This  correction  can  be  quite  large, 
an  order  of  magnitude  reduction  in  the  fatigue  life  at  low  frequencies,  compared 
to  high  frequency  results,  has  been  observed.  A  correction  should  also  be  made 
for  the  asymmetries  which  develop  in  thermal  cycling,  where  a  much  higher 
stress  is  reached  at  low  temperatures  than  at  high  temperatures  [5]. 

D.  Fatigue  Life  Criteria 

Assuming  one  has  calculated  A />  and  from  it  A7  and  from  it  Nn  one  is  now 
ready  to  define  the  fatigue  life.  The  question  now  arises,  however,  as  to  what 
constitutes  failure.  Failure  could  be  defined  as  the  number  of  cycles  N  i  to  grow 
a  crack  under  the  chip  carrier,  the  number  of  cycles  (N i  +  N2)  to  grow  it  into 
the  fillet  or  the  number  of  cycles  (N i  +  N2  +  N3)  to  grow  it  completely 
through  the  fillet.  Using  Nx  as  a  design  life  is  very  conservative,  probably  too 
conservative  since  if  the  fillet  is  uncracked  there  will  be  little  change  in  the  joint 
electrical  properties  [14].  Using  N\  +  N2  +  iV3  as  a  criteria  is  non  conserva¬ 
tive.  This  will  predict  the  average  behavior  but  it  has  been  shown  [13]  that  long 
before  the  average  joint  cracks  all  the  way  through  the  fillet  the  worst  joint  will 
have  already  failed  and  this,  not  the  average  joint  behavior,  marks  the  failure  of 
the  device. 

The  best  criteria  for  failure  is  a  resistance  increase  of  0.1%  in  the  first  joint 
to  show  such  an  increase.  Above  this  resistance  increase,  cracks  are  clearly 
observed  through  the  joint  fillet.  This  behavior  correlates  with  a  little  more  than 
a  50%  drop  in  load,  measured  in  the  joint  fatigue  tests.  This  in  turn  correlates 
with  N !  +  N2,  i.e.,  the  number  of  cycles  for  an  average  joint  to  have  a  crack 
growing  through  the  transition  region  into  the  fillet.  This  is  somewhat  conserva¬ 
tive  because  it  generally  takes  a  little  over  a  66%  drop  in  load  before  the  worst 
joint  reaches  a  0.1%  resistance  increase  and  a  visible  crack  through  the  fillet 
[13]. 
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In  summary,  the  suggested  design  approach  is  to 

A.  Calculate  the  thermal  displacement,  A p,  using  an  exact  solution  if  necessary. 

B.  Calculate  the  joint  strain  distribution  A-?,  from  A p  using  FEA. 

C.  Calculate  Nj  from  Ay,-  using  the  LCF  data,  correcting  for  frequency  and  asym¬ 
metry  effects. 

D.  Calculate  Nf  from  Nj,  correcting  for  the  joint  size  effect  described  here.  Define 
Nf  as  TNj,  Nf  =  N\  +  N 2  is  suggested,  where  Nx  +  N2  are  the  number  of 
cycles  to  grow  a  crack  into,  but  not  through  the  fillet  of  the  average  joint. 
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